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1. Introduction 
We consider the following problem. Let 

(1) P(t) = x n - a^x"- 1 + ■■■ + (-l) n a n (t) 

be a polynomial with all roots real, smoothly parameterized by t near in M. Can 
we find n smooth functions Xi(t), . . . , x n (t) of the parameter t defined near 0, which 
are the roots of P(t) for each tl We can reduce the problem to a\ = 0, replacing 
the variable x by with the variable y = x — ai(t)/n. We will say that the curve (1) 
is smoothly solvable near t — if such smooth roots Xi{t) exist. 
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We describe an algorithm which in the smooth and in the holomorphic case 
sometimes allows to solve this problem. The main results are: If all roots are real 
then they can always be chosen differentiable, but in general not C 1 for degree 
n > 3; and in degree 2 they can be chosen twice differentiable but in general not 
C 2 . If they are arranged in increasing order, they depend continuously on the 
coefficients of the polynomial, and if moreover no two of them meet of infinite order 
in the parameter, then they can be chosen smoothly. We also apply these results to 
obtain a smooth 1-parameter perturbation theorem for selfadjoint operators with 
compact resolvent under the condition, that no pair of eigenvalues meets of infinite 
order. 

We thank C. Fcffcrman who found a mistake in a first version of 2.4, to M. 
and Th. Hoffmann- Ostenhof for their interest and hints, and to Jerry Kazdan for 
arranging [15]. 

2. Choosing differentiable square and cubic roots 

2.1. Proposition. The case n = 2. Let P(t)(x) — x 2 — /(t) for a function f 
defined on an open interval, such that f{t) > for all t. 

If f is smooth and is nowhere flat of infinite order, then smooth solutions x exist. 

If f is C 2 then C 1 -solutions exist. 

If f is C 4 then twice differentiable solutions exist. 

Proof. Suppose that / is smooth. If /(to) > then we have obvious local smooth 
solutions ±y/ f(t). If /(to) = we have to find a smooth function x such that 
/ = x 2 , a smooth square root of /. If / is not flat at to then the first nonzero 
derivative at to has even order 2m and is positive, and f(t) = (t — t ) 2m f2 m (t), 

where f2m(t) ■= Jq ^(2m-i)\ f^ 2m \^o + r (t — to))dr gives a smooth function and 
hm(to) = j2^f (2m) (to) > 0. Then x(t) := (t - t ) m ^/hJf) is a local smooth 
solution. One can piece together these local solutions, changing sign at all points 
where the first non-vanishing derivative of / is of order 2m with m odd. These 
points are discrete. 

Let us consider now a function / > of class C 2 . We claim that then x 2 = f{t) 
admits a C 1 -root x(t), globally in t. We consider a fixed to- If /(to) > then there 
is locally even a C 2 -solution x±(t) = ±y/f(t). If /(to) = then /(t) = {t-t ) 2 h(t) 
where h > is continuous and C 2 off to with h(to) = \f"{t<s)- If M*o) > then 
x±(t) = ±(t - to)v/Mt) is C 2 off t , and 

x > ±{to) = Um X±(t |-; ±fa) = hm ±VM*) = ±vm = ±Jy(to). 

t—>to X — to £— *to v 

If h(t ) = then we choose x±(to) = 0, and any choice of the roots is then 
differentiable at to with derivative 0, by the same calculation. 

One can piece together these local solutions: At zeros t of / where f"(t) > 
we have to pass through 0, but where f"(t) = the choice of the root does not 
matter. The set {t : f(t) = f"(t) = 0} is closed, so its complement is a union of 
open intervals. Choose a point in each of these intervals where /(t) > and start 
there with the positive root, changing signs at points where /(t) = ^ /"(t): these 
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points do not accumulate in the intervals. Then we get a differentiable choice of 
a root x(t) on each of this open intervals which extends to a global differentiable 
root which is on {t : f(t) = f"(t) = 0}. 
Note that we have 

/'(*) 



At) 



m if/(*)>o 

±77W2 if /(*) = 



In points to with /(to) > the solution x is C 2 ; locally around points to with 
/(to) = and /"(to) > the root x is C 1 since for t ^ t near t wc have /(t) > 
and /'(t) ^ 0, so by l'Hospital we get 

,. // n2 ,. f'W 2 ,. 2/'(t)/"(t) /"(t ) „ N2 
lim x'(t) 2 = lim 44tt = lim 7777 n = "^T^ = x (*o) , 
t-t t-tc 4/(t) i-t 4/'(t) 2 v ; 

and since the choice of signs was coherent, x' is continuous at to; if /"(to) = then 
x'{to) = and x'(t) — > for t — > t for both expressions, by lemma 2.2 below. Thus 
x is C 1 . 

If moreover / > is C 4 , then the solution x from above may be modified 
to be twice differentiable. Near points to with /(to) > any continuous root 
t i ► x±(t) = ±y/ f(t) is even C 4 . Near points to with /(to) = /'(*o) = we have 
/(t) = (t — t ) 2 h(t) where h > is C 2 . We may choose a C 1 -root z with z 2 — h by 
the arguments above, and then x(t) := (t — t )z(t) is twice differentiable at to since 
we have 

x'(t) - x'(to) _ z(t) + (t - to)z'(t) - z(t ) 



t — to t — to 



If /(*o) = f"(to) = f^ 4 \to) = then any C 1 choice of the roots is twice differen- 
tiable at t , in particular x(t) = |t — to|z(t). 

Now we can piece together this solutions similarly as above. Let y be a global 
C 1 -root of /, chosen as above changing sign only at points t with f(t) = < /"(0). 
We put x(t) — e(t)y(t), where e(t) e {±1} will be chosen later. The set {t : /(t) = 
/"(t) = /( 4 )(t) = 0} has a countable union of open intervals as complements. In 
each of these intervals choose a point to with /(to) > 0, near which y is C 4 . Now 
let e(to) = 1, and let e change sign exactly at points with f(t) = f"(t) = but 
f^\t) > 0. These points do not accumulate inside the interval. Then x is twice 
differentiable. □ 

2.2. Lemma. Let f > be a C 2 -function with /(to) = 0, then for all t £ K we 

have 

(1) f'(t) 2 < 2/(t) max{/"(t + r(t - to)) : < r < 2}. 
Proof. If /(t) = then f'(t) = so (1) holds. We use the Taylor formula 

(2) f(t + s) = f(t) + f'(t)s + [ (1 - r)f"(t + rs) dr s 2 
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In particular we get (replacing t by to and then to + s by t) 

(3) f(t) = + 0+ f (1 - r)f"(t + r(t - to)) dr (t - t ) 2 

Jo 

< {t ~ 2 * 0) max{/"(t + r(t - to)) : < r < 2} 
Now in (2) we replace s by — es (where e = sign(/'(t))) to obtain 

(4) < f(t - es) = f(t) - \f'(t)\s + [ (1 - r)f"(t - ers) dr s 2 

Jo 

Let us assume t > to and then choose (using (3)) 



/ 2f(t) 

S{t) V max{/"(to + r(t-to)) : < r < 2} " * ~ * ' 

Note that we may assume /(t) > 0, then s(t) is well defined and s(t) > 0. This 
choice of s in (4) gives 

\2 



\f(t)\ < ± (f{t) + ^ max{/"(t - ers(t)) : < r < 1}^ 

(/(*) + ^ max{/"(i - r(t - ««,)) : -1 < r < 1}) 



= y/2 f(t) max{/"(t + r(t - to)) : < r < 2} 



s(t) 

which proves (1) for t > to. Since the assertion is symmetric it then holds for all 
t. □ 

2.3. Examples. If / > is only C 1 , then there may not exist a differentiable root 
of x 2 = f(t), as the following example shows: x 2 = f(t) := t 2 sin 2 (logi) is C 1 , but 
±tsin(logf) is not differentiable at 0. 

If / > is twice differentiable there may not exist a C 1 -root: x 2 = f(t) = 
i 4 sin 2 (j) is twice differentiable, but ±t 2 sin(^) is differentiable, but not C 1 . 

If / > is only C 3 , then there may not exist a twice differentiable root of 
x 2 = f(t), as the following example shows: x 2 = f(t) := t 4 sin 2 (log t) is C 3 , but 
±t 2 sin(log£) is only C 1 and not twice differentiable. 

2.4. Example. If f(t) > is smooth but flat at 0, in general our problem has 
no C 2 -root as the following example shows, which is an application of the general 
curve lemma 4.2.15 in [5]: Let h : M. — > [0, 1] be smooth with h(t) = 1 for t > and 
h(t) = for t < -1. Then the function 

°° /2n 1 \ 

f(t) : = h n(t - tn). ( ^r(* " tn? + ^ J , where 

n— 1 ^ ' 

M*): = ^(» 9 (r^fr+*))^(» a (-^T-*)) and 
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is > and is smooth: the sum consists of at most one summand for each t, and 
the derivatives of the summands converge uniformly to 0: Note that h n (t) = 1 for 
|*| < „. 2 Ui and h n (t) = for \t\ > n 2 i +i + 4y hence h n (t - t n ) ^ only for 
r n <t < r n+1 , where r n := YJl=\ (pr + ^+r)- Let c «( s ) : = f£ s2 + > and 
Hi := sup{|/t«(f)| : t G R}. Then 

n 2 sup{|(^. C „)( fe )(t)| :teR} = n 2 sup{\(h n -c n )^(t)\ : 1*1 < + ^} 

< n 2 £ (*Wsup{|#-0 (t) | : |t| < + -L } 

i=0 v ' 

^ (E (!) n2i+2 ^) SU P{I C ^ Wl : 1*1 ^ 2 > 3< k} 

and since c„ is rapidly decreasing in C°°(R, R) (i.e. {p(n) c n : n e N} is bounded in 
C°° (R, R) for each polynomial p) the right side of the inequality above is bounded 
with respect to n G N and hence the series ^2 n h n ( —t n )c n ( —t n ) converges 
uniformly in each derivative, and thus represents an clement of / e C°°(R, R). 
Moreover we have 

1 2T7 

f(t n ) = ^, f(tn) = 0, /"(«„) = —. 

Let us assume that /(i) = g(*) 2 for i near sup„ i„ < oo, where g is twice differen- 
tiablc. Then 

/' = 2.93' 

/" = 2. 9 .g" + 2(.g') 2 

2/r = 4 5 v' + (/') 2 

2/(tn)/"(tn) = ^{t n fg"{t n ) + /'(t„) 2 

thus g"(t n ) = ±2n, so g cannot be C 2 , and 5' cannot satisfy a local Lipschitz 
condition near limt„. Another similar example can be found in 7.4 below. 

According to [15], some results of this section are contained in Frank Warners 
dissertation (around 1963, unpublished): Non-negative smooth functions have C 1 
square roots whose second derivatives exist everywhere. If all zeros arc of finite 
order there are smooth square roots. However, there are examples not possessing 
a C 2 square root. Here is one: 

f(t) = sin 2 (l/i)e _1/ ' + e- 2/t for t > 0, f(t) = for t < 0. 

This is a sum of two non-negative C°° functions each of which has a C°° square 
root. But the second derivative of the square root of / is not continuous at the 
origin. 

In [6] Glaeser proved that a non-negative C 2 -function on an open subset of R" 
which vanishes of second order has a C 1 positive square root. A smooth function 
/ > is constructed which is flat at such that the positive square root is not C 2 . 
In [4] Dieudonne gave shorter proofs of Glaeser 's results. 
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2.5. Example. The case n > 3. We will construct a polynomial 

P(t) = x 3 + a 2 (t)x - a 3 (t) 

with smooth coefficients a 2 , a 3 with all roots real which does not admit C 1 -roots. 
Multiplying with other polynomials one then gets polynomials of all orders n > 3 
which do not admit C 1 -roots. 

Suppose that P admits C 1 -roots x\, x 2 , x 3 . Then we have 

= xi + x 2 + x 3 
a 2 = xix 2 + X2X3 + x 3 xi 
a 3 = xix 2 x 3 

= ±1 + x 2 + x 3 
a 2 = xix 2 + xix 2 + x 2 x 3 + x 2 x 3 + x 3 xi + x 3 i\ 
a 3 = X1X2X3 + X1X2X3 + xix 2 x 3 

We solve the linear system formed by the last three equations and get 

a 3 - a 2 xi 



xi = 
x 2 = 
x 3 = 



(x 2 - x 1 )(x 3 - x-l) 
a 3 - a 2 x 2 

(x 3 - x 2 )(xi - x 2 ) 
a 3 - a 2 x 3 



(xi - x 3 )(x 2 - x 3 )' 

We consider the continuous function 

... af, + a 2 a 3 a 2 — a 2 a 3 

b 3 := x x x 2 x 3 = ,3,0-2 • 

4a 2 + 27 a 3 

For smooth functions / and e with e 2 < 1 we let 

u : = -I2a 2 := f 
v : = 108a 3 := sf. 

then all roots of P are real since a 2 < and 432(4a|+27a§) = v 2 -u 3 = f 6 (e 2 -l) < 
0. We get then 



116646 3 



Now we choose 



4i> 3 - 27uii 2 v + 27ii 3 v 



f 3 e 3 + 9f 2 fe 2 e + 27 f f 2 e(e 2 - 1) + 27f 3 e{e 2 - 3) 
e 2 - 1 



OO 

/(*): = 2 M* -*«)■(£(*-*«)). 

n=l 

1 



n=l 
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where h n and t n are as in the beginning of 2.4. Then f(t n ) — 0, f(t n ) = jk, and 
£(t n ) — grr, hence 



1086 3 (fn) 



f(t n ) 3 e(t n )(e(t n ) 2 -3) 
e{tn) 2 - 1 



oo 



is unbounded on the convergent sequence t n . So the roots cannot be chosen locally 
Lipschitz, thus not C x . 

3. Choosing local roots of real polynomials smoothly 

3.1. Preliminaries. We recall some known facts on polynomials with real coeffi- 
cients. Let 

P(x) = x n - axx"- 1 + • • • + (-l)"a„ 

be a polynomial with real coefficients oi, . . . , a n and roots x\, . . . , x n € C. It 
is known that ai = ai(xi, . . . , x n ), where Oi (i — 1, ...,n) are the elementary 
symmetric functions in n variables: 



ai(xi, 



E 



l<ji<---<ji<n 



Denote by Si the Newton polynomials Y^j=i x \ which are related to the elementary 
symmetric function by 



(I) Sfc — Sfc-ifi + Sfc-202 H h (-1) sier fe _i + (-l) K fccr fc = (fc < n) 

The corresponding mappings are related by a polynomial diffeomorphism ip n , given 
by (1): 



a 

s % 
s r 



= (o-i, 

= (si, 
= f 1 o o-" 



Sn) ■ 



Note that the Jacobian (the determinant of the derivative) of s n is n\ times the 
Vandermonde determinant: det(ds n (x)) — nl Yii>j( x i — x j) ~- n ' Van(x), and even 
the derivative itself d(s n )(x) equals the Vandermonde matrix up to factors i in the 
i-th row. We also have det(d(ip n )(x)) = (-l)"(" +3 )/ 2 n! = (-I)"^- 1 )/ 2 ™!, and 
consequently det(da n (x)) — Y\ i:> j(xj — Xi). We consider the so-called Bezoutiant 



B := 



/ s 
si 



si 

S2 



Sn-1 \ 



V S n ... S2n-2 I 

Let Bk be the minor formed by the first k rows and columns of B. From 



B k (x) 



t 1 

Xl 



1 

X2 



I \ f I Xl 

1 x 2 



Jc-1 



ct 1 

„fe-l 



\ 



\1 a;« 



„fc-i 
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it follows that 

(2) A fc (x) := det(B fe (a;)) = (x n -x l2 ) 2 . . . (x h -x in ) 2 . . . (x lk _ 1 ~x lk ) 2 , 

il<i2<---<ik 

since for n x /c-matrices A one has det{AA T ) = J2 il< ... <ik det(A ily ,,^ ik ) 2 , where 
Ai lt ...,i k is the minor of A with indicated rows. Since the are symmetric we 
have Afe = Afe o <j n for unique polynomials A& and similarly we shall use B. 

3.2. Theorem. (Sylvester's version of Sturm's theorem, see [14], [10] J The roots 
of P are all real if and only if the symmetric (n x n) matrix B(P) is positive 
semidefinite; then Afc(P) := Afc(ai, . . . , a n ) > for 1 < k < n. The rank of B 
equals the number of distinct roots of P and its signature equals the number of 
distinct real roots. 

3.3. Proposition. Let now P be a smooth curve of polynomials 

P(t)(x) =x n - a^ty 1 - 1 + ■■■ + (-l) n a n (t) 

with all roots real, and distinct for t = 0. Then P is smoothly solvable near 0. 

This is also true in the real analytic case and for higher dimensional parameters, 
and in the holomorphic case for complex roots. 

Proof. The derivative ^P(0)(x) does not vanish at any root, since they are dis- 
tinct. Thus by the implicit function theorem we have local smooth solutions x(t) 
of P(t,x) = P(t)(x) = 0. □ 

3.4. Splitting Lemma. Let P be a polynomial 

P (x) - x n - a x x n - x + ■■■ + (-l) n a n . 

If Po = Pi ■ P2 , where Pi and Pi are polynomials with no common root. Then for P 
near Pq we have P = Pi(P)P2(P) for real analytic mappings of monic polynomials 
P 1 ► Pi{P) and P 1 ► PiiP), defined for P near P , with the given initial values. 

Proof. Let the polynomial Po be represented as the product 

Po = P1.P2 = (x p - hx^ 1 + ■■■ + (-l) p b p )(x« - ci**- 1 + • • • + (-1)%), 

let Xi for i = 1, . . . ,n be the roots of Po, ordered in such a way that for i — 1, . . . ,p 
we get the roots of Pi, and for i = p + 1, . . . ,p + q = n we get those of Pi. Then 
(ai, . . . , a n ) — (j?> q {bi, . . . , b p , a, . . . , c q ) for a polynomial mapping <j) p ' q and we get 

a n = 4?' q o (a p x a g ), 
det (da n ) = det(d<p p > q (b,c))det(do- p )det(d<r q ). 

From 3.1 we conclude 

Y\ (xi - xj) = dct(d<p p ' q (b,c)) Y\ ( x i~ x j) I] i^i-Xj) 

l<i<j<n l<i<j<p p+l<i<j<n 

which in turn implies 

det(d^ q (b, c)) = Yl ( x i - x i)^ {) 

l<i<p<j<n 

so that (j) 9 ^ is a real analytic diffeomorphism near (6, c). □ 
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3.5. For a continuous function / defined near in R let the multiplicity or order 
of flatness m(f) at be the supremum of all integer p such that f(t) = t p g{t) 
near for a continuous function g. If f is C n and m(f) < n then f(t) = t m ^g(t) 
where now g is C n ~ m ^ and g(0) ^ 0. If / is a continuous function on the space 
of polynomials, then for a fixed continuous curve P of polynomials we will denote 
by m(f) the multiplicity at of t ^ f(P(t)). 

The splitting lemma 3.4 shows that for the problem of smooth solvability it is 
enough to assume that all roots of P(0) are equal. 

Proposition. Suppose that the smooth curve of polynomials 

P(t)(x) = x n + a 2 {t)x n - 2 -■■■ + (-l) n a n (t) 

is smoothly solvable with smooth roots t i— > Xi{t), and that all roots of P(0) are 
equal. Then for (k = 2, . . . , n) 

fn(Ak) > k(k — 1) min m(xi) 

\<i<n 

m(afc) > k min m{xi) 

l<i<n 

This result also holds in the real analytic case and in the smooth case. 

Proof. This follows from 3.1.(2) for Afc, and from a,k(t) — <7k(%i{t), ■ ■ ■ i x n {t)). □ 

3.6. Lemma. Let P be a polynomial of degree n with all roots real. If a\ = a 2 = 
then all roots of P are equal to zero. 

Proof. From 3.1.(1) we have £ x\ = s 2 (x) = af(x) - 2a 2 (x) =a\- 2a 2 = 0. □ 

3.7. Multiplicity lemma. Consider a smooth curve of polynomials 

P(t)(x) =x n + a 2 {t)x n - 2 -■■■ + {-l) n a n (t) 

with all roots real. Then, for integers r, the following conditions are equivalent: 

(1) m(afe) > kr for all 2 < k < n. 

(2) m(Afc) > k(k - l)r for all2<k<n. 

(3) m(a 2 ) > 2r. 

Proof. We only have to treat r > 0. 

(1) implies (2): From 3.1.(1) we have m(sfe) > rk, and from the definition of 
Afc = det(B fc ) we get (2). 

(2) implies (3) since A2 = — 2na 2 . 

(3) implies (1): From 02(0) = and lemma 3.6 it follows that all roots of the 
polynomial P(0) are equal to zero and, then, 03(0) = • • • = a„(0) = 0. There- 
fore, m(as) , ■ ■ ■ , m(a n ) > 1. Under these conditions, we have a 2 (t) = t 2r a 2t2r (t) 
and ak(t) = t mk ak,m k {t) for k = 3, . . . , n, where the are positive integers and 
fl 2,2r, 013,77131 • ■ • 7 dn,m n are smooth functions, and where we may assume that either 
TOfe = m(dfc) < 00 or TOfc > kr. 

Suppose now indirectly that for some k > 2 we have = m(afe) < kr. Then 
we put 

. . m 3 m„, 
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We consider the following continuous curve of polynomials for t > 0: 

P m (t)(x) := x n + a 2 . 2r {t)t 2r - 2m x n - 2 

- a 3 , m3 (t)t m3 - 3m x n - 3 + ■■■ + (-l) n a n , mn (t)t m "- nm . 

If xi, . . . , x n are the real roots of P(t) then t~ m xi, . . . , t~ m x n are the roots of P m (t), 
for t > 0. So for t > 0, P m {t) is a family of polynomials with all roots real. Since 
by theorem 3.2 the set of polynomials with all roots real is closed, P m (0) is also a 
polynomial with all roots real. 

By lemma 3.6 all roots of the polynomial P m (0) are equal to zero, and for 
those k with m,k — km we have Ofc iTOfe (0) = and, therefore, m(afc) > m,k, a 
contradiction. □ 

3.8. Algorithm. Consider a smooth curve of polynomials 

P{t){x) = x n - ai^x™- 1 + a 2 (t)x n - 2 + {-l) n a n (t) 

with all roots real. The algorithm has the following steps: 

(1) If all roots of P(0) are pairwise different, P is smoothly solvable for t near 
by 3.3. 

(2) If there are distinct roots at t — we put them into two subsets which splits 
P(t) = Pi(t).P 2 {t) by the splitting lemma 3.4. We then feed P l (t) (which 
have lower degree) into the algorithm. 

(3) All roots of -P(O) are equal. We first reduce P(t) to the case a\(t) — by 
replacing the variable x by y = x — ai(t)/n. Then all roots are equal to 
so m(a 2 ) > 0. 

(3a) If 771(02) is finite then it is even since A2 = —2na 2 > 0, m(a 2 ) = 2r and by 
the multiplicity lemma 3.7 Oj(t) = ai.i r (t)t' tr (i = 2, . . . , n) for smooth a,^. 
Consider the following smooth curve of polynomials 

P r (t)(x) =x n + a 2i2r (t)x n - 2 - a 3i3r (t)x n - 3 ■■■ + (-l) n a n , nr (t). 

If P r (t) is smoothly solvable and Xk(t) are its smooth roots, then Xk(t)t r are 
the roots of P(t) and the original curve P is smoothly solvable too. Since 
fl2,2m(0) 7^ 0, not all roots of P r (0) are equal and we may feed P r into step 
2 of the algorithm. 

(3b) If 771(02) is infinite and 02 = 0, then all roots are by 3.6 and thus the 
polynomial is solvable. 

(3c) But if 771(02) is infinite and 02 7^ 0, then by the multiplicity lemma 3.7 all 
m(aj) for 2 < i < n arc infinite. In this case we keep P(t) as factor of the 
original curve of polynomials with all coefficients infinitely flat at t — after 
forcing oi = 0. This means that all roots of P(t) meet of infinite order of 
flatness (see 3.5) at t = for any choice of the roots. This can be seen as 
follows: If x(t) is any root of P(t) then y(t) := x(t)/t r is a root of P r (t), 
hence by 4.1 bounded, so x(t) = t r ^ 1 .ty{t) and t 1— > ty(t) is continuous at 
t = 0. 
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This algorithm produces a splitting of the original polynomial 

P(t) =p(°°)(t)pW(t) 

where p(°°) has the property that each root meets another one of infinite order at 
t = 0, and where P( s \t) is smoothly solvable, and no two roots meet of infinite 
order at t — 0, if they are not equal. Any two choices of smooth roots of P^ differ 
by a permutation. 

The factor p(°°) may or may not be smoothly solvable. For a flat function / > 
consider: 

x 4 - (f(t) + t 2 )x 2 + t 2 f(t) = (x 2 - f(t)).(x - t)(x + t). 

Here the algorithm produces this factorization. For f(t) = g(t) 2 the polynomial is 
smoothly solvable. For the smooth function / from 2.4 it is not smoothly solvable. 

4. Choosing global roots of polynomials differentiably 
4.1. Lemma. For a polynomial 

P(x) =x n - oi(P)a; n - 1 + • • • + (-l)"a„(P) 

with all roots real, i.e. Afc(P) = Afe(oi, . . . , a n ) > for 1 < k < n, let 

xi(P) <x 2 (P) < ••• <z„(P) 

be the roots, increasingly ordered. 

Then all x% : a n (M. n ) — > R are continuous. 

Proof. We show first that x\ is continuous. Let Po € cr"(IR n ) be arbitrary. We have 
to show that for every e > there exists some 5 > such that for all \P — Pq\ < 6 
there is a root x(P) of P with x(P) < xi(P ) + e and for all roots x(P) of P we have 
x(P) > xi(Pq) — e. Without loss of generality we may assume that xi(Pq) = 0. 

We use induction on the degree n of P. By the splitting lemma 3.4 for the C°-case 
we may factorize P as Pi(P)-p2(P), where Pi(Po) has all roots equal to x\ = and 
^MPj) has all roots greater than and both polynomials have coefficients which 
depend real analytically on P. The degree of p2(P) is now smaller than n, so by 
induction the roots of p2(P) are continuous and thus larger than xi(P ) — e for P 
near P . 

Since was the smallest root of Po we have to show that for all e > there exists 
a S > such that for |P — P | < 5 any root x of Pi(P) satisfies \x\ < e. Suppose 
there is a root x with \x\ > e. Then we get as follows a contradiction, where n\ is 
the degree of Pi. From 

- xni = J2(-^ ka ^ xni ~ k 
fe=i 

we have 

ni n\ n± ^ 

e<\x\ = ^(-lfafctPi)^! < M p i)l ^ < Y. £ — £l ~ k = £ > 

fc=l k=\ k=\ 711 
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provided that n\\ak{Pi)\ < £ k , which is true for Pi near P , since ak(Po) = 0. Thus 
x\ is continuous. 

Now we factorize P = (x — x\{P)) ■ P 2 (P), where P 2 (P) has roots x 2 (P) < 
•■■ < x n (P). By Horner's algorithm (o„ = b n -\X\, a n -i = b n -i + b n - 2 Xi, . . . , 
a 2 = b 2 + biXi, a\ =bi+x\) the coefficients bk of P 2 (P) are again continuous and 
so we may proceed by induction on the degree of P. Thus the claim is proved. □ 

4.2. Theorem. Consider a smooth curve of polynomials 

P(t){x) =x n + a 2 (t)x n - 2 -■■■ + (-l) n a n {t) 

with all roots real, for t £ R. Let one of the two following equivalent conditions be 
satisfied: 

(1) If two of the increasingly ordered continuous roots meet of infinite order 
somewhere then they are equal everywhere. 

(2) Let k be maximal with the property that Afc(P) does not vanish identically 
for all t. Then Afc(P) vanishes nowhere of infinite order. 

Then the roots of P can be chosen smoothly, and any two choices differ by a 
permutation of the roots. 

Proof. The local situation. We claim that for any t , without loss to — 0, the 
following conditions are equivalent: 

(1) If two of the increasingly ordered continuous roots meet of infinite order at 
t = then their germs at t = are equal. 

(2) Let k be maximal with the property that the germ at t = of Afc(P) is not 
0. Then A fe (P) is not infinitely flat at t = 0. 

(3) The algorithm 3.8 never leads to step (3c). 

(3) ==> (1). Suppose indirectly that two nonequal of the increasingly ordered 
continuous roots meet of infinite order at t = 0. Then in each application of step (2) 
these two roots stay with the same factor. After any application of step (3a) these 
two roots lead to nonequal roots of the modified polynomial which still meet of 
infinite order at t = 0. They never end up in a facter leading to step (3b) or 
step (1). So they end up in a factor leading to step (3c). 

(1) (2). Let ari(i) < • • • < x n (t) be the continuous roots of P(t). From 3.1, (2) 
we have 

(4) A fc (P(t))= J2 (x tl -x l2 ) 2 ...(x n ~x l J 2 ...(x lk _ 1 -x lk ) 2 . 

ii<i2<---<ik 

The germ of A^(P) is not 0, so the germ of one summand is not 0. If Afe(P) were 
infinitely flat at t = 0, then each summand is infinitely flat, so there are two roots 
among the Xi which meet of infinite order, thus by assumption their germs are 
equal, so this summand vanishes. 

(2) =^> (3). Since the leading Afc(P) vanishes only of finite order at zero, P has 
exactly k different roots off 0. Suppose indirectly that the algorithm 3.8 leads to 
step (3c), then P = p(°°)p( s ) for a nontrivial polynomial p(°°). Let x x {t) < ■ ■ ■ < 
x p (t) be the roots of p(°°)(t) and x p+1 (t) < ■ ■ ■ < x n (t) those of P^~> . We know 



CHOOSING ROOTS OF POLYNOMIALS SMOOTHLY 



13 



that each x% meets some xj of infinite order and does not meet any xi of infinite 
order, for i,j<p<l. Let fc(°°) > 2 and be the number of generically different 
roots of p(°°) and p( s \ respectively. Then k = fc(°°) + k^ s \ and an inspection of 
the formula for Afc(P) above leads to the fact that it must vanish of infinite order 
at 0, since the only non- vanishing summands involve exactly fc(°°) many generically 
different roots from p(°°). 

Let 2/i (i) < • • • < Uk(t) be the From 3.1, thctag2 we 

The global situation. From the first part of the proof we see that the algorithm 
3.8 allows to choose the roots smoothly in a neighborhood of each point tel, and 
that any two choices differ by a (constant) permutation of the roots. Thus we may 
glue the local solutions to a global solution. □ 

4.3. Theorem. Consider a curve of polynomials 

P(t)(x) =x n - a^x 71 ' 1 + ■■■ + {-l) n a n {t), t e R, 

with all roots real, where all ai are of class C n . Then there is a differentiable curve 
x = (x\, . . . , x n ) : R — > M™ whose coefficients parameterize the roots. 

That this result cannot be improved to C 2 -roots is shown already in 2.4, and 
not to C 1 for n > 3 is shown in 2.5. 

Proof. First we note that the multiplicity lemma 3.7 remains true in the C"-case 

for r = 1 in the following sense, with the same proof: 

If a\ = then the following two conditions are equivalent 

(1) <Zfc(i) = t k ak,k(t) for a continuous function ak,k, for 2 < k < n. 

(2) a%{t) — t 2 a2 l 2(t) for a continuous function 02.2- 

In order to prove the theorem itself we follow one step of the algorithm. First 
we replace x by x + ^a\(t), or assume without loss that a\ = 0. Then we choose a 
fixed t, say t = 0. 

If 02(0) = then it vanishes of second order at 0: if it vanishes only of first 
order then A 2 (P(i)) = — 2na 2 (t) would change sign at t = 0, contrary to the 
assumption that all roots of P(t) are real, by 3.2. Thus 02(f) = t 2 a2^{t) 1 so by 
the variant of the multiplicity lemma 3.7 described above we have ak(t) = t k ak,k(t) 
for continuous functions ak,k, for 2 < fc < n. We consider the following continuous 
curve of polynomials 

Pi(t)(x) =x n + a 2a {t)x n - 2 - a 3 . 3 (t)x n - 3 ■■■ + (-l) n a n , n (t). 

with continuous roots z\(t) < ■■■ < z n (t), by 4.1. Then Xk(t) = Zk(t)t are dif- 
ferentiable at 0, and are all roots of P, but note that Xk(t) — yk(t) for t > 0, 
but Xk(t) — y n -k(t) for t < 0, where yi(t) < ••• < y n {t) arc the ordered roots of 
P(t). This gives us one choice of differentiable roots near t — 0. Any choice is then 
given by this choice and applying afterwards any permutation of the set {1, . . . , n} 
keeping invariant the function k Zk(0). 

If 02(0) 7^ then by the splitting lemma 3.4 for the C™-case we may factor 
P(t) = P\{t) . . . Pk(t) where the Pi(t) have again C™-coefficients and where each 
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Pj(0) has all roots equal to Cj, and where the Cj are distinct. By the arguments 
above the roots of each Pj can be arranged differentiably, thus P has differentiable 
roots y k (t). 

But note that we have to apply a permutation on one side of to the original 
roots, in the following case: Two roots Xk and xi meet at zero with Xk(t) — xi(t) = 
tcki(t) with Cfcz(O) ^ (we say that they meet slowly). We may apply to this choice 
an arbitrary permutation of any two roots which meet with Cki(0) — (i.e. at least 
of second order), and we get thus any differentiable choice near t = 0. 

Now we show that we choose the roots differentiable on the whole domain R. 
We start with with the ordered continuous roots y\ (t) < • • • < y n (t) . Then we put 

Xk{t) = Va(t){k){t) 

where the permutation <j(t) is given by 

a{t) = (1, 2) £1 > 2 W ... (1, n) ei -»W (2, 3) £2 - 3(t) . . . (n - 1, nf-^ 

and where £ij{t) <G {0, 1} will be specified as follows: On the closed set Sij of all t 
where yi (t) and yj (t) meet of order at least 2 any choice is good. The complement 
of Si_j is an at most countable union of open intervals, and in each interval we 
choose a point where we put Ei.j = 0. Going right (and left) from this point we 
change Eij in each point where yi and yj meet slowly. These points accumulate 
only in Sij. □ 

5. The real analytic case 
5.1. Theorem. Let P be a real analytic curve of polynomials 

P{t){x) =x n - ai^z"- 1 + • • • + (-l) n a n (t), t G R, 

with all roots real. 

Then P is real analytically solvable, globally on M. All solutions differ by per- 
mutations. 

By a real analytic curve of polynomials we mean that all ai(t) are real analytic 
in t (but see also [8]), and real analytically solvable means that we may find Xi{t) 
for i = 1, . . . ,n which are real analytic in t and are roots of P(t) for all t. The 
local existence part of this theorem is due to Rellich [11], Hilfssatz 2, his proof uses 
Puiseux-expansions. Our proof is different and more elementary. 

Proof. We first show that P is locally real analytically solvable, near each point 
to E M. It suffices to consider to = 0. Using the transformation in the introduction 
we first assume that a\(t) = for all t. We use induction on the degree n. If n = 1 
the theorem holds. For n>lwc consider several cases: 

The case 02(0) 7^ 0. Here not all roots of P(0) are equal and zero, so by the 
splitting lemma 3.4 we may factor P(t) = Pi(t).P 2 (t) for real analytic curves of 
polynomials of positive degree, which have both all roots real, and we have reduced 
the problem to lower degree. 
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The case a 2 (0) = 0. If a 2 (t) = for all t, then by 3.6 all roots of P(t) are 0, 
and we are done. Otherwise 1 < m(a 2 ) < oo for the multiplicity of a 2 at 0, and 
by 3.6 all roots of P(0) arc 0. If m(a 2 ) > is odd, then A 2 (P)(t) = -2na 2 {t) 
changes sign at t = 0, so by 3.2 not all roots of P(t) are real for t on one side of 0. 
This contradicts the assumption, so m(a 2 ) = 2r is even. Then by the multiplicity 
lemma 3.7 we have ai(t) = cns r {t)t lr (i = 2, . . . ,n) for real analytic a ijir , and we 
may consider the following real analytic curve of polynomials 

P r (t)(x) = x n + a 2 . 2r {t)x n - 2 - a 3 . 3r (t)x n - 3 ■■■ + {-l) n a n , nr (t), 

with all roots real. If P r (t) is real analytically solvable and Xk(t) are its real analytic 
roots, then Xk{t)t r are the roots of P(t) and the original curve P is real analytically 
solvable too. Now a 2 , 2r (0) ^ and we are done by the case above. 

Claim. Let x = (xi, . . . , x n ) : I — > R™ be a real analytic curve of roots of P on 
an open interval I C R. Then any real analytic curve of roots of P on I is of the 
form a o x for some permutation a. 

Let y : I — > R™ be another real analytic curve of roots of P. Let tk — > to be a con- 
vergent sequence of distinct points in /. Then y(tk) = ctk{x{tk)) — {x ak i, ■ ■ ■ , x a k n) 
for permutations a^. By choosing a subsequence we may assume that all are 
the same permutation a. But then the real analytic curves y and a o x coincide on 
a converging sequence, so they coincide on / and the claim follows. 

Now from the local smooth solvability above and the uniqueness of smooth so- 
lutions up to permutations we can glue a global smooth solution on the whole of 
R. □ 

5.2. Remarks and examples. The uniqueness statement of theorem 5.1 is wrong 
in the smooth case, as is shown by the following example: x 2 = f(t) 2 where / is 
smooth. In each point t where / is infinitely flat one can change sign in the solution 
x(t) = ±f(t). No sign change can be absorbed in a permutation (constant in t). If 
there are infinitely many points of flatness for / we get uncountably many smooth 
solutions. 

Theorem 5.1 reminds of the curve lifting property of covering mappings. But 
unfortunately one cannot lift real analytic homotopies, as the following example 
shows. This example also shows that polynomials which arc real analytically pa- 
rameterized by higher dimensional variables are not real analytically solvable. 

Consider the 2-parameter family x 2 = t\ + 1 2 . The two continuous solutions are 
x(t) = ±\t\, but for none of them t\ ^> x(t\, 0) is diffcrentiable at 0. 

There remains the question whether for a real analytic submanifold of the space 
of polynomials with all roots real one can choose the roots real analytically along 
this manifold. The following example shows that this is not the case: 

Consider 

P(h,t 2 )(x) = {x 2 - (t 2 + 4)) (x - (h - oi)) {x - (t 2 - a 2 )), 

which is not real analytically solvable, see above. For ai ^ a 2 the coefficients 
describe a real analytic embedding for (ti,t 2 ) near 0. 
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6. Choosing roots of complex polynomials 

6.1. In this section we consider the problem of finding smooth curves of complex 
roots for smooth curves t P(t) of polynomials 

P(t)(z) =z n - a^tjz"- 1 + ■■■ + (-l) n a n (t) 

with complex valued coefficients ai(t), . . . , a n (t). We shall also discuss the real 
analytic and holomorphic cases. The definition of the Bezoutiant B, its principal 
minors Afe, and formula 3.1.(2) are still valid. Note that now there are no restric- 
tions on the coefficients. In this section the parameter may be real and P(t) may 
be smooth or real analytic in t, or the parameter t may be complex and P(t) may 
be holomorphic in t. 

6.2. The case n = 2. As in the real case the problem reduces to the following 
one: Let / be a smooth complex valued function, defined near in M, such that 
/(0) = 0. We look for a smooth function g : (K, 0) — ► C such that / = g 2 . If m(/) is 
finite and even, we have f(t) = t m ^h{t) with h(0) ^ 0, and g(t) := t m ^l 2 ^/t[fy 
is a local solution. If m(f) is finite and odd there is no solution g, also not in the 
real analytic and holomorphic cases. If f(t) is flat at t = 0, then one has no definite 
answer, and the example 2.4 is still not smoothly solvable. 

6.3. The general case. Proposition 3.3 and the splitting lemma 3.4 are true in 
the complex case. Proposition 3.5 is true also because it follows from 3.1.(2). Of 
course lemma 3.6 is not true now and the multiplicity lemma 3.7 only partially 
holds: 

6.4. Multiplicity Lemma. Consider the smooth (real analytic, holomorphic) 
curve of complex polynomials 

P(t)(z) = z n + a 2 {t)z n - 2 -■■■ + (-l) n a n (t). 

Then, for integers r, the following conditions are equivalent: 

(1) m(afe) > kr for all 2 < k < n. 

(2) m(Afc) > k(k - l)r for all 2 < k < n. 

Proof. (2) implies (1): Since A2 = na 2 we have £2(0) = — 2a2(0) = 0. From 
As(O) = — S3(0) 2 we then get s 3 (0) = 0, and so on we obtain s 4 (0) = • • • = s„(0) = 
0. Then by 3.1.(1) aj(0) = for i — 3, . . . , n. The rest of the proof coincides with 
the one of the multiplicity lemma 3.7. □ 

6.5. Algorithm. Consider a smooth (real analytic, holomorphic) curve of poly- 
nomials 

P(t)(z) =z n - ax^z"- 1 + a 2 (t)z n - 2 -■■■ + (-l)"a„(i) 

with complex coefficients. The algorithm has the following steps: 

(1) If all roots of P(0) are pairwise different, P is smoothly (real analytically, 
holomorphically) solvable for t near by 3.3. 

(2) If there are distinct roots at t — we put them into two subsets which 
factors P(t) = Pi(t).P 2 (t) by the splitting lemma 3.4. We then feed P t (t) 
(which have lower degree) into the algorithm. 
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(3) All roots of P(0) are equal. We first reduce P(t) to the case a\(t) — by 
replacing the variable x by y = x — ai(t)/n. Then all roots are equal to 
so a,i(0) = for all i. 

If there does not exist an integer r > with m{ai) > ir for 2 < i < n, then 
by 3.5 the polynomial is not smoothly (real analytically, holomorphically) 
solvable, by proposition 3.5: We store the polynomial as an output of the 
procedure, as a factor of P^> below. 

If there exists an integer r > with m(ai) > ir for 2 < i < n, let 
a,i(t) = ai^ r {t)t %r (i = 2, ...,n) for smooth (real analytic, holomorphic) 
ai^i r . Consider the following smooth (real analytic, holomorphic) curve of 
polynomials 



If P r (t) is smoothly (real analytically, holomorphically) solvable and Xk(t) 
are its smooth (real analytic, holomorphic) roots, then Xk(t)t r are the roots 
of P{t) and the original curve P is smoothly (real analytically, holomorphi- 
cally) solvable too. 

If for one coefficient we have m{ai) = ir then P r (0) has a coefficient 
which does not vanish, so not all roots of P r (0) are equal, and we may feed 
P r into step (2). 

If all coefficients of P r (0) are zero, we feed P r again into step (3). 
In the smooth case all m(a,i) may be infinite; In this case we store the 
polynomial as a factor of p(°°) below. 

In the holomorphic and real analytic cases the algorithm provides a splitting of 
the polynomial P(t) — P^ [t)P^ sS> (t) into holomorphic and real analytic curves, 
where P( s \t) is solvable, and where P( n \t) is not solvable. But it may contain 
solvable roots, as is seen by simple examples. 

In the smooth case the algorithm provides a splitting near t = 



the same property as above. 

6.6. Remarks. If P(t) is a polynomials whose coefficients are meromorphic func- 
tions of a complex variable t, there is a well developed theory of the roots of 
P(t)(x) = as multi- valued meromorphic functions, given by Puiseux or Laurent- 
Puiseux series. But it is difficult to extract holomorphic information out of it, and 
the algorithm above complements this theory. See for example Theorem 3 on page 
370 (Anhang, §5) of [1]. The question of choosing roots continuously has been 
treated in [2]: one finds sufficient conditions for it. 

7. Choosing eigenvalues and eigenvectors of matrices smoothly 

In this section we consider the following situation: Let A(t) = (Aij(t)) be a 
smooth (real analytic, holomorphic) curve of real (complex) (n x n)-matrices or 



P r (t)(x) =x n + a 2 ^ r {t)x n - 2 - a 3 , 3r (t)x 



,71 — 3 



■ ■ ■ + (-l) n On,nr(t)- 
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operators, depending on a real (complex) parameter t near 0. What can we say 
about the eigenvalues and eigenfunctions of A(t)7 

Let us first recall some known results. These have some difficulty with the 
interpretation of the eigenprojections and the eigennilpotents at branch points of 
the eigenvalues, see [7], II, 1.11. 

7.1. Result. ([7], II, 1.8) Let C3(h A(t) be a holomorphic curve. Then all 
eigenvalues, all eigenprojections and all eigennilpotents are holomorphic with at 
most algebraic singularities at discrete points. 

7.2. Result. ([11], Satz 1) Let t i— > A(t) be a real analytic curve of hermitian com- 
plex matrices. Let X be a k-fold eigenvalue of A{Q) with k orthonormal eigenvectors 
Vi, and suppose that there is no other eigenvalue of A near X. Then there are k real 
analytic eigenvalues Xi(t) through X, and k orthonormal real analytic eigenvectors 
through the Vi, for t near 0. 

The condition that A(t) is hermitian cannot be omitted. Consider the following 
example of real semisimple (not normal) matrices 



where at t = we do not get a base of eigenvectors. 

7.3. Result. (Rellich [13], see also Kato [7], II, 6.8) Let A(t) be a C 1 -curve of 
symmetric matrices. Then the eigenvalues can be chosen C 1 in t, on the whole 
parameter interval. 

For an extension of this result to Hilbert space, under stronger assumptions, 
see 7.8, whose proof will need 7.3. This result is best possible for the degree of 
continuous differentiability, as is shown by the following example. 

7.4 Example. Consider the symmetric matrix 



The characteristic polynomial of A(t) is A 2 — (a(t) 2 + b(t) 2 ). We shall specify the 
entries a and b as smooth functions in such a way, that a(t) 2 + b(t) 2 does not admit 
a C 2 -square root. 

Assume that a{t) 2 + b(t) 2 = c(t) 2 for a C 2 -function c. Then we may compute as 
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follows: 



2 2,72 

c = a + b 



cc = 



aa' + bb' 



(c') 2 + cc" = {a'f + aa" + {b'f + bb" 



l{{a'f + aa" + {b'f + bb"-{c'f) 

l -{(a'f + aa" + m 2 + bb"- {aa '\ hh ' )2 
c\ c 

(aV - ba') 2 + ttV + b 3 b" + ab 2 a" + a?bb" 
c 3 



Bye 2 



6 2 we have 



< 1, 



b 3 



< 1, 



< 



V3' 



a 2 6 



< 



V3' 



So for C 2 -functions a, b, and continuous c all these terms are bounded. We will 
now construct smooth a and b such that 



A2 



(gfr' - fa 7 ) 

C 3 



(a 2 + 6 2 ) 3 



is unbounded near t = 0. This contradicts that c is C 2 . For this we choose a and 6 
similar to the function / in example 2.4 with the same t n and h n : 

a(t) : = £ Mi - *n)- (^(t - t n ) + ±) , 



n=l 



K*) : = SM* -*»)■( ^r(*- *»»))■ 



n=l 



2'^ 
2n . 



Then o(t„) = ±, b(t n ) - 0, |c(t„)| - ±, and fe'(t„) - |?. 

7.5. Result. (Rellich [12], see also Kato [7], VII, 3.9) Let A(t) be a real analytic 
curve of unbounded self-adjoint operators in a Hilbert space with common domain 
of definition and with compact resolvent. 

Then the eigenvalues and the eigenvectors can be chosen real analytically in t, 
on the whole parameter domain. 

7.6. Theorem. Let A(t) = (Aij(t)) be a smooth curve of complex hermitian 
(n x n)-matrices, depending on a real parameter iel, acting on a hermitian space 
V = C n , such that no two of the continuous eigenvalues meet of infinite order at 
any tell/ they are not equal for all t. 

Then all the eigenvalues and all the eigenvectors can be chosen smoothly in t, on 
the whole parameter domain K. 

The last condition permits that some eigenvalues agree for all t — we speak of 
higher 'generic multiplicity' in this situation. 
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Proof. The proof will use an algorithm. 

Note first that by 4.2 the characteristic polynomial 

(1) P(A(t))(\) = det(A(t) - XI) 

= A" - a^A"- 1 + a 2 (i)A"- 2 - • • • + (-l) n a n (t) 

n 

= ^Trace(AM(t))A™-* 

i=0 

is smoothly solvable, with smooth roots Ai(i), . . . A n (i), on the whole parameter 
interval. 

Case 1: distinct eigenvalues. If A(0) has some eigenvalues distinct, then one can 
reorder them in such a way that for io = < 1 < i\ < it < ■ ■ ■ < ik < n = iu+i we 
have 

Ai(0) = ■ ■ ■ = K (0) < A n+1 (0) = • • • = A i2 (0) < • • • < A ife+ i(0) = • • • = A n (0) 
For t near we still have 

Ai(t),...,A 4l (t) < X ll+1 (t),...,X l2 (t) < ... < X ik+1 (t),...,X n (t) 
Now for j = 1, . . . , k + 1 consider the subspaces 

V} j) = {v€V:(A(t)-X i (t))v = 0} 

i=ij_i + l 

Then each V t ^ runs through a smooth vector subbundle of the trivial bundle 
(— e,e) x V — > (—£,£), which admits a smooth framing e^_ 1+ i(t), . . . , e^. (t). We 

have V = V t U) for each t. 

In order to prove this statement note that 

V® = kcr((A(t) - A^._ 1+1 (i)) o . . . o (A(t) X t] (t))) 

so is the kernel of a smooth vector bundle homomorphism B{t) of constant rank 
(even of constant dimension of the kernel), and thus is a smooth vector subbundle. 
This together with a smooth frame field can be shown as follows: Choose a basis of 
V such that A(0) is diagonal. Then by the elimination procedure one can construct 
a basis for the kernel of B(Q). For t near 0, the elimination procedure (with the 
same choices) gives then a basis of the kernel of B(t); the elements of this basis are 
then smooth in t, for t near 0. 

From the last result it follows that it suffices to find smooth eigenvectors in each 
subbundle V^) separately, expanded in the smooth frame field. But in this frame 
field the vector subbundle looks again like a constant vector space. So feed each of 
this parts (A restricted to as matrix with respect to the frame field) into case 
2 below. 
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Case 2: All eigenvalues at are equal. So suppose that A(t) : V — > V is hermitian 
with all eigenvalues at t = equal to ai ^ , see (1). 

Eigenvectors of A(t) arc also eigenvectors of A(t) — ^3-^-1, so we may replace 

A(t) by A(t) - ^1 and assume that for the characteristic polynomial (1) we have 
«i = 0, or assume without loss that Aj(0) = for all i and so A(0) = 0. 
If A(t) = for all t we choose the eigenvectors constant. 

Otherwise let Aij(t) = tA^\t). From (1) we see that the characteristic polyno- 
mial of the hermitian matrix A^(t) is P\{t) in the notation of 3.8, thus m(aj) > i 
for 2 < i < n (which follows from 3.5 also). 

The eigenvalues of A^> (t) are the roots of P\ (t) , which may be chosen in a smooth 
way, since they again satisfy the condition of theorem 4.2. Note that eigenvectors 
of A^ are also eigenvectors of A. If the eigenvalues are still all equal, we apply 
the same procedure again, until they are not all equal: we arrive at this situation 
by the assumption of the theorem. Then we apply case 1. 

This algorithm shows that one may choose the eigenvectors Xi(t) of A(t) in a 
smooth way, locally in t. It remains to extend this to the whole parameter interval. 

If some eigenvalues coincide locally, then on the whole of R, by the assumption. 
The corresponding eigenspaces then form a smooth vector bundle over R, by case 1, 
since those eigenvalues, which meet in isolated points arc different after application 
of case 2. 

So we we get V — where each is a smooth sub vector bundles of 

V x R, whose dimension is the generic multiplicity of the corresponding smooth 
eigenvalue function. It suffices to find global orthonormal smooth frames for each 
of these vector bundles; this exists since the vector bundle is smoothly trivial, by 
using parallel transport with respect to a smooth Hermitian connection. □ 

7.7. Example, (see [11], §2) That the last result cannot be improved is shown by 
the following example which rotates a lot: 

(II)- M*)==(-^), M*)=±e-*, 
A{t) : = (*+(*),*_(*)) ( A+ W A _° (t) ) (s+W.M*))- 1 
= e-^( C ° S l S HV 

V sin t - cos t J 

Here t A(t) and t A±(t) are smooth, whereas the eigenvectors cannot be 
chosen continuously. 

7.8. Theorem. Let t i— > A(t) be a smooth curve of unbounded self-adjoint op- 
erators in a Hilbert space with common domain of definition and with compact 
resolvent. Then the eigenvalues of A(t) may be arranged in such a way that each 
eigenvalue is C 1 . 

Suppose moreover that no two of the continuously chosen eigenvalues meet of 
infinite order at any t € R if they are not equal. Then the eigenvalues and the 
eigenvectors can be chosen smoothly in t, on the whole parameter domain. 
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Remarks. That A(t) is a smooth curve of unbouded operators means the follow- 
ing: There is a dense subspace V of the Hilbert space H such that V is the domain 
of definition of each A(t), and such that A{t)* = A(t) with the same domains V, 
where the adjoint operator A(t)* is defined by (A(t)u,v) = (u,A(t)*v) for all v for 
which the left hand side is bounded as function in u e H. Moreover we require that 
t i ► (A(t)u,v) is smooth for each u e V and v e H. This implies that t i— ► A(t)u 
is smooth R -» i? for each u e V by [9], 2.3 or [5], 2.6.2. 

The first part of the proof will show that t A(t) smooth implies that the 
resolvent (A(t) — z)^ 1 is smooth in t and z jointly, and only this is used later in 
the proof. 

It is well known and in the proof we will show that if for some (t, z) the resolvent 
(A(t) — z)^ 1 is compact then for all £ e R and z in the resolvent set of A(t). 

Proof. For each t consider the norm ||u|| 2 := ||w|| 2 + ||^4(£)u|| 2 on V. Since A(t) = 
A(t)* is closed, (V, \\ \\t) is again a Hilbert space with inner product (u,v) t := 
(u, v) + (A(t)u, A(t)v) . All these norms are equivalent since (V, | ||t + || \\ s ) — * 
(V, || ||t) is continuous and bijective, so an isomorphism by the open mapping 
theorem. Then t i— > (u,v) t is smooth for fixed u, v G V, and by the multilinear 
uniform boundedness principle ([9], 5.17 or [5], 3.7.4 + 4.1.19) the mapping t 
( , } t is smooth into the space of bounded bilinear forms on (V", || || s ) for each 
fixed s. By the exponential law ([9], 3.12 or [5], 1.4.3) (t,u) i— > ||w|| 2 is smooth 
from R x (V, || || s ) — » K for each fixed s. Thus all Hilbert norms || ||t are locally 
uniformly equivalent, since {||u|| t : \t\ < K, \\u\\ s < 1} is bounded by Lk, s in R, so 
\\u\\t < iif,s||u||s for all \t\ < K. Let us now equip V with one of the equivalent 
Hilbert norms, say || ||o- Then each A(t) is a globally defined operator V — > H 
with closed graph and is thus bounded, and by using again the ( multilinear uniform 
boundedness theorem as above we see that 1 1— » A(t) is smooth R — > L(V, iJ). 

If for some (i, z) e 1 x C the bounded operator A{t) — z : V — > i? is invcrtiblc, 
then this is true locally and (t, z) i— > (A(t) — z) _1 : i? — > V is smooth since inversion 
is smooth on Banach spaces. 

Since each A(t) is hermitian the global resolvent set {(t, z) e R x C : (A(t) — z) : 
V — > H is invertible} is open, contains R x (C \ R), and hence is connected. 

Moreover (A(t) — z)^ 1 : H — > H is a compact operator for some (equivalently 
any) (t, z) if and only if the inclusion i : V — > i? is compact, since i = (A(t) — z) _1 o 



Let us fix a parameter s. We choose a simple smooth curve 7 in the resolvent 
set of A(s) for fixed s. 

(1) Claim For i near s, there are C 1 -functions t Aj(t) : 1 < i < N which 
parametrize all eigenvalues (repeated according to their multiplicity) of A(t) 
in the interior of 7. If no two of the generically different eigenvalues meet 
of infinite order they can be chosen smoothly. 

By replacing A(s) by A(s)—zq if necessary we may assume that is not an eigenvalue 
of A(s). Since the global resolvent set is open, no eigenvalue of A{t) lies on 7 or 
equals 0, for t near s. Since 



(A(t) — z) : V —> H —> H. 
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is a smooth curve of projections (on the direct sum of all eigenspaces corresponding 
to eigenvalues in the interior of 7) with finite dimensional ranges, the ranks (i.e. 
dimension of the ranges) must be constant: it is easy to see that the (finite) rank 
cannot fall locally, and it cannot increase, since the distance in L(H,H) of P(i) to 
the subset of operators of rank < N — rank(F(s)) is continuous in t and is either 
or 1. So for t near s, there are equally many eigenvalues in the interior, and we 
may call them ^(t) : 1 < i < N (repeated with multiplicity). Let us denote by 
£i(t) ■ 1 < i < N a corresponding system of eigenvectors of A(t). Then by the 
residue theorem we have 

N 

mityeiWieiit), ) = -tt- z P (Mt) - z)- 1 dz 
i=i zm J ~t 

which is smooth in t near s, as a curve of operators in L(H, H) of rank N, since 
is not an eigenvalue. 

(2) Claim. Let t 1— > T(t) £ L(H, H) be a smooth curve of operators of rank N in 
Hilbert space such that T(0)T(0)(H) = T(0)(H). Then t h-> Trace(T(t)) is 
smooth near (note that this implies T smooth into the space of nuclear op- 
erators, since all bounded linear functionals are of the form A 1— ► Tracc(AF), 
by [9], 2.3 or 2.14.(4). 

Let F := T(0)(H). Then T(t) = (T^t), T 2 (t)) : H -> F F 1 - and the image of 
T(t) is the space 

T(t)(H) = {(T^),^ (*)(*)) :xeH} 

= {(T 1 (t)(x),T 2 (t)(x)) : x e F} for t near 

= {{y, S(t)(y)) :yeF}, where S(t) := T 2 (t) o (Ti^F)" 1 . 

Note that S(t) : F — > F 1 - is smooth in t by finite dimensional inversion for T\(t)\F : 
F -> F. Now 

™)=- e ((_^ gg>j£)(^ ;)) 

- Trace ( ( T M\ F T ^\ F± 

-M( r r r,( r)u„ :))— 

= Tracc ^(t)|F + (T^f -^(i) TxWIF- 1 - 

- Tracc(Ti(t)|F+ (Ti(t)|F J -)S(t) : F^ f), 

which is visibly smooth since F is finite dimensional. 
From claim (2) we now may conclude that 

N 

^2^{t) p = -—Trace J z p {A{t) - z)' 1 dz 



1 
S(t) 1 
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is smooth for t near s. 

Thus the Newton polynomial mapping s N (fii(t), . . . , /Uiv(i)) i s smooth, so also 
the elementary symmetric polynomial <r N (fii(t), . . . , /Uiv(i)) is smooth, and thus 
{Hi{t) : 1 < i < N} is the set of roots of a polynomial with smooth coefficients. 
By theorem 4.3 there is an arrangement of these roots such that they become 
diffcrcntiable. If no two of the generically different ones meet of infinite order, by 
theorem 4.2 there is even a smooth arrangement. 

To see that in the general case they are even C 1 note that the images of the 
projections P(t, 7) of constant rank for t near s describe the fibers of a smooth 
vector bundle. The restriction of A{t) to this bundle, viewed in a smooth framing, 
becomes a smooth curve of symmetric matrices, for which by Rellich's result 7.3 
the eigenvalues can be chosen C x . This finishes the proof of claim (1). 

(3) Claim. Let t Xi(t) be a differentiable eigenvalue of A(t), defined on some 
interval. Then 

|A 4 (i 1 )-A,(i 2 )|<(l + |A,(i 2 )|)( e a l tl - t2 l-l) 

holds for a continuous positive function a = a(ii,i 2 ) which is independent 
of the choice of the eigenvalue. 

For fixed t near s take all roots Xj which meet A^ at t, order them diffcrcntiably near 
t, and consider the projector P(t, 7) onto the joint cigenspaces for only those roots 
(where 7 is a simple smooth curve containing only Xi(t) in its interior, of all the 
eigenvalues at t). Then the image of u <—> P(u, 7), for u near t, describes a smooth 
finite dimensional vector subbundle of M x H, since its rank is constant. For each u 
choose an othonormal system of eigenvectors Vj(u) of A(u) corresponding to these 
Xj(u). They form a (not necessarily continuous) framing of this bundle. For any 
sequence tk — > t there is a subsequence such that each vj(tk) — > Wj(t) where Wj(t) 
is again an orthonormal system of eigenvectors of A(t) for the eigenspace of Aj(i). 
Now consider 

A{t) - X t (t) A{t k ) - A(t) \(t k )-X t (t) 

— 7 7 — ViW H -, -, Vi(tk) Vi(t k ) = 0, 

tk — 1 ife — t t k — t 

take the inner product of this with Wi(t), note that then the first summand vanishes, 
and let tk — > t to obtain 

A-(i) = (A'(t)wi(t) 7 Wi(t)) for an eigenvector Wi(t) of A(t) with eigenvalue Xi(t). 

This implies, where Vt = (V, || \\t), 

\Xm < \\A'(t)\\ L{VtM) \\ Wl (t)\\v t \\Mt)\\H 

= \\A'{t)\\Hv t ,H)yJ\\wi(t)\\% + \\mMt)\\ 2 H 
= \\A'(t)\\ L{Vt , H) y/1 + Ai(i) 2 < a + a\Xi(t)\, 

for a constant a which is valid for a compact interval of t's since i 1 — s- 1 1 \\^ is 
smooth on V. By GronwalPs lemma (see e.g. [3], (10.5.1.3)) this implies claim (3). 
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By the following arguments we can conclude that all eigenvalues may be num- 
bered as Xi(t) for i in N or Z in such a way that they are differentiable (by which 
we mean C 1 , or C°° under the stronger assumption) in t e M. Note first that by 
claim (3) no eigenvalue can go off to infinity in finite time since it may increase at 
most exponentially. Let us first number all eigenvalues of ^4(0) increasingly. 

We claim that for one eigenvalue (say Ao(0)) there exists a differentiable extension 
to all of R; namely the set of all £ e M with a differentiable extension of Ao on the 
segment from to t is open and closed. Open follows from claim (1). If this intervall 
does not reach infinity, from claim (3) it follows that (t, Ao(t)) has an accumulation 
point (s,x) at the the end s. Clearly x is an eigenvalue of A(s), and by claim (f) 
the eigenvalues passing through (s, x) can be arranged differentiably, and thus \o(t) 
converges to x and can be extended differentiably beyond s. 

By the same argument we can extend iteratively all eigenvalues differentiably to 
all (el: if it meets an already chosen one, the proof of 4.3 shows that we may 
pass through it coherently 

Now we start to choose the eigenvectors smoothly, under the stronger assump- 
tion. Let us consider again eigenvalues {Xi(t) : 1 < i < N} contained in the interior 
of a smooth curve 7 for t in an open interval /. Then V t := P(t, 7) (if) is the fiber 
of a smooth vector bundle of dimension N over I. We choose a smooth framing 
of this bundle, and use then the proof of theorem 7.6 to choose smooth sub vector 
bundles whose fibers over t are the eigenspaces of the eigenvalues with their generic 
multiplicity. By the same arguments as in 7.6 we then get global vector sub bundles 
with fibers the eigenspaces of the eigenvalues with their generic multiplicity, and 
finally smooth eigenvectors for all eigenvalues. □ 
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